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If cold dark matter is present at the galactic center, as in current models of the dark halo, it is 
accreted by the central black hole into a dense spike. Particle dark matter then annihilates strongly 
inside the spike, making it a compact source of photons, electrons, positrons, protons, antiprotons, 
and neutrinos. The spike luminosity depends on the density profile of the inner halo: halos with 
finite cores have unnoticeable spikes, while halos with inner cusps may have spikes so bright that 
the absence of a detected neutrino signal from the galactic center already places interesting upper 
limits on the density slope of the inner halo. Future neutrino telescopes observing the galactic center 
could probe the inner structure of the dark halo, or indirectly find the nature of dark matter. 
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The evidence is mounting for a massive black hole at 
the galactic center. Ghez et al. jjj have confirmed and 
sharpened the Keplerian behavior of the star velocity dis- 
persion in the inner 0.1 pc of the galaxy found by Eckart 
and Genzel || . These groups estimate the mass of the 
black hole to be M = 2.6 ± 0.2 x 10 6 M Q . 

If cold dark matter is present at the galactic center, 
as in current models of the dark halo, it is redistributed 
by the black hole into a cusp. We call it the central 
'spike,' to avoid confusion with the inner halo cusp fa- 
vored by present N-body simulations of galaxy forma- 
tion H . If cold dark matter contains neutral elementary 
particles that can annihilate with each other, like the 
supersymmetric neutralino, the annihilation rate in the 
spike is strongly increased as it depends on the square of 
the matter density. The steep spike profile, with index 
> 3/2, then implies that most of the annihilations take 
place at the inner radius of the spike, determined either 
by self-annihilation or by capture into the black hole. 

Of the annihilation end-products, neutrinos escape the 
spike and propagate to us undisturbed. Current limits on 
the neutrino emission from the galactic center place up- 
per limits on the slope of the inner halo. Future neutrino 
telescopes may improve on these limits or bring informa- 
tion on the nature of dark matter. 



I. ADIABATIC SPIKE AROUND THE CENTRAL 
BLACK HOLE 

We find the dark matter density profile in the re- 
gion where the black hole dominates the gravitational 
potential. From the data in 0,0], this is the region 
t < Rm — 0.2 pc. Other masses (the central star cluster, 
for example) also influence the dark matter distribution, 



but since they make the gravitational potential deeper, 
their effect is to increase the central dark matter density 
and the annihilation signals. 

We work under the assumption that the growth of the 
black hole is adiabatic. This assumption is supported by 
the collisionless behavior of particle dark matter. We can 
find the final density after black hole formation from the 
final phase-space distribution f'(E',L') as 
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We have neglected the contribution from unbound orbits 
(E' > 0). The lower limit of integration L' c , and the sec- 
ond factor in E' m , are introduced to eliminate the parti- 
cles with L < 2cRs which are captured by the black hole. 
(Rq = 2GM/c 2 is the Schwarzschild radius.) We relate 
f'(E', L') to the initial phase-space distribution f(E, L) 
through the relations, valid under adiabatic conditions, 
/'(£', £/) = f(E, L), V = L, I'(E', V) = I{E, L), where 
the last two equations are the conservation of the angular 
momentum L and of the radial action I(E, L). 

The density slope in the spike depends not only on the 
slope of the inner halo but also on the behavior of the 
initial phase-space density f(E,L) as E approaches the 



potential at the center (f)(0) m. If f(E,L) approaches 
a constant, as in models with finite cores, the spike 
slope is 7 S p = 3/2. If f(E,L) diverges as [E - 4)(0)}- fJ , 
as in models with an inner cusp, the spike slope is 
7s P > 3/2. Models with finite cores include ||: the 
non-singular and the modified isothermal sphere, the 
Henon isochrone model, the Plummer model, the King 
models, the modified Hubble profile, the Evans power- 
law models with R c ^ 0, and the Zhao (a, (3, 7) models, 
p ~ r -7(i + r i/«)-/3a ; with 7 = and 1/(2q) = integer. 
Models with an inner cusp include 0]: the models of 
Jaffe, Hernquist, Navarro-Frenk-White, the 7/7/ models 
of Dehnen and Tremaine et al., and the other Zhao 
models. 

As an example of models with finite cores, we consider 
the isothermal sphere. It has f(E,L) = pQ(2ira 2 )~ 3 / 2 
cxp(—E/a 2 ). Close to the black hole, we have E <C ct 2 , 
and f(E,L) ~ po(2na 2 )^ 3 / 2 , a constant. Then from 
eq. (||) we easily find 
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valid for r <C Rm — 0.2 pc. The last factor comes from 
the capture of particles by the black hole: the density 
vanishes for r < 4i?s- 

As an example of models with an inner cusp, we 
consider a single power law density profile, p(r) = 
Po{ r / r o)~ 1 ', with < 7 < 2. Its phase-space distribu- 
tion function is 
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with 13 = (6- 7 )/[2(2- 7 )] and O = 4irGr 2 p / '[(3- 7 )(2 - 
7)]. To find f'(E',L'), we need to solve I'[E',L') = 
I(E, L) for E as a function of E' . In the field of a point- 
like mass, we have I'(E', L') = 2tt [-L 1 + GM/^TE 1 ] . 
In the field of the power law profile, whose potential is 
proportional to r 2-7 , the action integral cannot be per- 
formed exactly. We have found an approximation good 
to better than 8% over all of phase space for < 7 < 2: 
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where A = [2/(4 - 7 )]V(2-7) [(2 - 7 )/(4 - 7 )] 1 /2 an d b 
7r(2- 7 )/B( 5 i 7 ,§). Expressing,/ 
and integrating eq. (QJ), we obtain 



7r(2 — j)/B(jl—, |). Expressing £ as a function of E' 
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with p R = po (R sp /r a ) 7 , 7 sp = (9 - 27)/ (4 - 7), and 
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For < 7 < 2, the density 



slope in the spike, 7 sp , varies only between 2.25 and 2.5. 
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FIG. 1. Examples of spike density profiles. 

While the exponent 7 sp can also be obtained by scaling 
arguments B , the normalization a 7 and the factor g 1 (r) 
accounting for capture must be obtained numerically. We 
find that g 7 (r) ~ (1 — 4i?s/r) 3 over our range of 7, and 
that a 7 ~ 0.2937 4 / 9 for 7 < 1, and is a 7 = 0.00733, 
0.120, 0.140, 0.142, 0.135, 0.122, 0.103, 0.0818, 0.0177 at 
7 = 0.05, 0.2, 0.4, . . . , 1.4, 2. The density falls rapidly to 
zero at r < 9.55Rg, vanishing for r < 4i?g. 

Annihilations in the inner regions of the spike set a 
maximal dark matter density p COIC — m/avtbht where 
tbh is the age of the black hole, conservatively 10 10 yr, 
m is the mass of the dark matter particle, and av is its 
annihilation cross section times relative velocity (notice 
that for non-relativistic particles av is independent of v). 
Using dp/ dt — —avp 2 /m, the final spike profile is 
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which has a core of radius R COTC — R sp (pn/ Pcore) ' ' sp - In 
the particle models we consider, not more than the initial 
amount of dark matter within 300 pc is annihilated. 

Examples of spike density profiles are shown in Fig. 1. 

To conclude this section, wc derive a conservative esti- 
mate of the dark matter density near the galactic center. 
Assume that the halo density is constant on concentric 
ellipsoids. Then the halo contribution Vh(r) to the ro- 
tation speed at distance r fixes the halo mass within r. 
Assume further that the halo density decreases monoton- 
ically with distance. Since at large radii it decreases at 
least as fast as r -2 , and at small radii only as r -7 with 
7 < 2, the density profile becomes steeper with distance. 
So to continue the r -7 dependence to all radii keeping 
the same halo mass interior to r as given by the rotation 
speed, we must decrease the density normalization. In 
this way we obtain an underestimate of the density near 
the center. Letting po — po(D/ro) _7 , we have 
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where we have taken Vh(D) = 90kms _1 at the Sun dis- 
tance D = 8.5 kpc, as obtained after subtracting a (some- 
what overestimated) luminous matter contribution of 180 
km s _1 to the circular speed of 220 ± 20kms~ 1 11. 



II. SIGNALS FROM NEUTRALINO 
ANNIHILATIONS 

Our analysis applies in general to a self-annihilating 
dark matter particle. To make it concrete, we examine 
a case of supersymmetric dark matter, the lightest neu- 
tralino. The minimal supersymmetric standard model 
privides a well-defined calculational framework, but con- 
tains at least 106 yet- unmeasured parameters H. Most 
of them control details of the squark and slepton sectors, 
and can safely be disregarded in dark matter studies. 
So we restrict the number of parameters to 7, follow- 
ing Bergstrom and Gondolo |g] . Out of the database of 
points in parameter space built in refs. |p|-|lO||, we use the 
35121 points in which the neutralino is a good cold dark 
matter candidate, in the sense that its relic density sat- 
isfies 0.025 < fl x h 2 < 1. The upper limit comes from the 
age of the Universe, the lower one from requiring that 
neutralinos are a major fraction of galactic dark halos. 

Gravitational interactions bring the cold neutralinos 
into our galactic halo and into the central spike, where 
neutralino pairs can annihilate and produce photons, 
electrons, positrons, protons, antiprotons, and neutrinos. 
While most products are subject to absorption and/or 
diffusion, the neutrinos escape the spike and propagate 
to us undisturbed. We focus on the neutrinos and post- 
pone the study of other signals. 

The expected neutrino flux from neutralino annihila- 
tions in the direction of the galactic center can be divided 
into two components: emission from the halo along the 
line of sight and emission from the central spike, 
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The halo flux from neutralino annihilations between us 
and the galactic center can be estimated assuming that 
a single power law profile p(r) — puir/D) -7 extends out 
to the Sun position r = D. The integrated neutrino flux 
within an angle O of the galactic center is 
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where m is the neutralino mass, Y v is the number of 
neutrinos produced per annihilation, either differential 
or integrated in energy, av is the neutralino-neutralino 
annihilation cross section times relative velocity, and 
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Here O is in radians and O m i„ 
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We evaluate the neutrino yield Y v and the neutralino 
annihilation cross sections a using the DarkSUSY pack- 
age [0], which incorporates Pythia simulations of the 
v continuum jl2| and the annihilation cross sections 
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FIG. 2. Enhancement of annihilation signals from the 
galactic center. 
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FIG. 3. For a Moore et al. halo profile, flux of neu- 
trino-induced muons in a neutrino telescope from neutralino 
dark matter annihilations in the direction of the galactic cen- 
ter, with (upper panel) and without (lower panel) the central 
spike. The horizontal line is the current upper limit. 

To the halo flux we need to add the contribution from 
the spike around the black hole at the galactic center. 
For an isothermal distribution we find 
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the factor of 25 serving to match the exact integration. 
This flux is a factor of ~ 10~ 9 smaller than the flux from 
dark matter annihilations between us and the galactic 
center, and so the addition of the spike does not modify 
the signal. The same conclusion is reached in general for 
halo models with finite cores. 

A strong enhancement results instead for power law 
profiles. We find 
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where R sp is given after eq. (M) with po replaced by 
Pd and ro by D. We fix i?i n so as to match the inte- 



gration of the numerically-calculated density profile in- 
cluding capture and annihilation: we find that R ln = 

I try 

f.5 [(20i?s) 2 + -R? re] gives a good approximation to 
the flux (within 6% for our values of 7). 

Contrary to the case of finite cores, for cusped halos 
there is a huge increase in flux from the galactic center 
when the spike is included, typically 5 orders of magni- 
tude or more, unless the inner halo slope 7 is very small 
(see Fig. 2, where = 1?5). The enhancement is no- 
table, for example, for the profile of Moore et al. which 
has 7 = 1.4 (see Fig. 3): including the spike around the 
black hole dramatically changes the prospects of observ- 
ing a neutrino flux. 

The neutrino flux from the spike increases with the 
inner halo slope 7. Imposing that it does not exceed 
the observed upper limit of ~ 2 x 10 4 muons(> lGcV) 
km -2 yr _1 [[L4| leads to an upper bound on 7. There is a 
separate upper bound for each model. They are plotted 
in Fig. 4a. (Plotted values of 7 max > 2 are unphysical 
extrapolations but are shown for completeness.) Present 
bounds are of the order of 7 max ~ 1.5. 

Future neutrino telescopes situated in the Northern 
hemisphere may improve on this bound or find a signal. 
For example, with a muon energy threshold of 25 GeV, 
the neutrino flux from the spike after imposing the cur- 
rent constraints could still be over 3 orders of magnitude 
above the atmospheric background (Fig. 5), allowing to 
probe 7 as low as 0.05 (Fig. 4b). 
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FIG. 4. Maximum inner slope 7 of the dark matter halo 
compatible with the upper limit on the neutrino emission from 
the galactic center, (a) Current limit at 1 GeV; (b) future 
reach at 25 GeV. 

In conclusion, we have shown that if the galactic dark 
halo is cusped, as favored in recent N-body simulations of 
galaxy formation, a bright dark matter spike would form 
around the black hole at the galactic center. A search of 
a neutrino signal from the spike could either set upper 
bounds on the density slope of the inner halo or clarify 
the nature of dark matter. 

This research has been supported in part by grants 
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FIG. 5. Maximal flux of neutrino-induced muons in a neu- 
trino telescope from neutralino annihilations at the galactic 
center, after imposing the current constraints on the neutrino 
emission. 
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